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The theoretical description of a Raman amplifier based on the vector model of 
randomly birefringent fibers is proposed and applied to the characterization of 
Raman polarizers. The Raman polarizer is a special type of Raman amplifier with 
the property of producing a highly repolarized beam when fed by relatively weak 
and unpolarized light. © 2010 Optical Society of America 

OCIS codes: 230.5440; 060.4370; 230.1150; 230.4320 

Polarization-dependent gain (PDG), an intrinsic characteristic of optical fiber-based Raman am- 
plifiers, is generally considered an unwanted feature for telecom-related applications. Very recently 
such opinion about the role of PDG was reversed, as the quest for higher transmission capacities 
brings to the forefront the need for polarization multiplexing protocols and polarization-controlling 
devices. Indeed, Martinelli et. al. demonstrated in Ref. [1] such a device, called Raman polarizer, 
which selectively amplifies only one polarization mode of the input beam, and thereby yields only 
this mode at the output, independently of the input state of polarization (SOP) of the signal beam. 
The development of a simple, yet rigorous as well as computer-friendly theory of Raman polarizers 
along with the scheme for their characterization is thus the purpose of this Letter. 

Telecom fibers are randomly birefringent fibers. Representative examples of vector theories of 
Raman amplifiers developed for telecom fibers can be found in Refs. [2,3]. The analytic theory 
of Ref. [2] is limited by the condition that the beat length Lb is smaller than the birefringence 
correlation length L c , and therefore their validity is questionable when applied to Raman polarizers, 
which as we shall see require the opposite inequality Lb S> L c . The full-scale numerical approach in 
Ref. [3] accurately models a randomly birefringent fiber as consisting of fiber spans with randomly 
distributed values and orientations of the birefringence. Typically, thousands of such realizations are 
required for getting an accurate statistics. Hence the required computer time is three to four orders 
of magnitude longer than for the numerical modeling involved in the theory presented below. In 
addition to the much faster performance, our theory is formulated in terms of a set of deterministic 
differential equations, and as such allows for a simple physical interpretation. 
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Starting with the equations of motion formulated by Lin and Agrawal in Ref. [2] we extend the 
one-beam model of the stochastic fiber proposed by Wai and Menyuk in Ref. [4] to two beams 
interacting not only via Kerr, but also via Raman effect. Detailed derivations can be found in 
Ref. [5], while here we only provide the final equation formulated for the Stokes vector = 
(s[ s \ 63 , 1S3 ) of the signal beam: 

7(w s ) x JW(*)SW +SW x j x ( z )s^ 

+e s g (sfthmSW + S ( Q s) J R (z) S&) . (1) 

The components of the Stokes vector are written in terms of the two polarization components V^i 
and V S 2 of the slowly varying signal field in the appropriate reference frame, as Sf = V S \V* 2 + 
V*iV S 2, S% = i(V* 1 V S 2 — V s i V* 2 ), S3 = \V s i\ 2 — I V^2 1 2 - Similar equations and definitions (with 
labels p and s interchanged) hold for the pump beam. 7(w s ) is the Kerr coefficient of the fiber 
at frequency u s of the signal beam; go is the Raman gain coefficient; f3'(u s ) is the inverse group 
velocity of the signal beam; a s is the attenuation coefficient; e s = 1; e p = —u> p /u> s . Matrices in the 
Eq. (1) are all diagonal with elements Jr = diag( Jr\, Jr2, Jr.3), J x = diag(Jxi, Jx2, Jxz), J s = 
diag(Jsi, J52, ^53)- Here J B1 = (Re(tt? 4 - u? )), Jr 2 = -(Re(it? 4 + uf )), J R3 = ~{\uu\ 2 - \u 10 \ 2 ), 
Jxi = |(Re(uf + u\ 3 -u\- uf 4 )), J X 2 = §(Re(tt? + u? 4 - u% - u? 3 )), J X 3 = fd^P + |«u| 2 - 
\uvs\ 2 - \uio\ 2 }, Jsi = 5(Re(ujj)), Js2 = -§(Re(u§)), J53 = § [3(«|) - l], and also Jr = (|^io| 2 + 
|iti 4 | 2 ). The three groups of coefficients {(u±), (u 2 ), (u 2 )}, {(Re 2 (ii 4 )), (Re 2 («5)), (Re 2 (tig))}, and 
{(Im 2 (u 4 )), (Im 2 (u5)), (Im 2 (ti6))} obey equations 

8 t Gi = -2L- 1 (G 1 -G 2 ), 

d z G 2 = 2L^\G 1 - G 2 ) - 4A/3K)G 4 , 

8 Z G 3 = 4A/3(w s )G 4 , 

Z G 4 = -L~ l G 4 + 2A(3(u s )(G 2 - G 3 ) , 

when we associate them with {Gi, G2, G3} respectively. Initial con- 
ditions are respectively (1, 0, 0), (0, 1, 0), and (0, 0, 1). In turn, the 
rest four groups of coefficients {(Re 2 (u7)), (Re 2 (iis)), (Re 2 (itg), (Re 2 (uio))}, 
{<Im 2 ( U7 )>, <Im 2 ( U8 )>, (Im 2 (u 9 ), (Im 2 ( Ul0 ))}, {(Re 2 (u n )>, (Re 2 (u 12 )), (Re 2 (u 13 ), (Re 2 (u 14 ))}, 
and {(Im 2 (un)), (Im^-u^)), (Im (^13), (Im 2 (ui 4 ))}, can be found from equations 

d z G x = -2L~ 1 (G 1 - G 2 ) + 2A„G 5 , 

d z G 2 = 2L-\G l - G 2 ) - 2A+G 6 , 

8 Z G 3 = 2A+G 6 , 

8 Z G A = -2A_G 5 , 

8 Z G 5 = A_(G 4 -G 1 )-L- 1 G 5 , 

9 Z G 6 = A + (G 2 -G 3 )-L- 1 G 6 , 
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when we associate them with {G±, G2, G3, G4}, with initial conditions as (1, 0, 0, 0), (0, 0, 0, 1), 
(0, 1, 0, 0), and (0, 0, 1, 0), respectively. Here, A± = A/3(u) p ) ± A/3(w s ), where A(3(u) s ) [A/3(u p )] is 
the magnitude of the birefringence at frequency oj s (cj p ). The power of the signal beam defined as 



i2 /•_\2 /' i a2\ ^ 2 



S ( Q s) = ( S[ s) + S { 2 s) + S ( 3 sy 1 obeys the equation 



(d t + P'{u s )d t ) S { S) = -a s S { s) + go[JmS^S^ 

+J R1 S[ S) S^ + J R2 S { 2 s) si p) + J R3 si s) si p) ) . (2) 



Eqs. (1) and (2) for the signal (and pump) fields are the key finding of our study. These equations 
are valid for a wide range of parameters and regimes, for undepleted as well as with a depleted 
pump. The only limitation is that the total length of the fiber L and/or the nonlinear length 
Lnl = [i(. u; s )Sq' ] _1 be longer than the correlation length L c . Eqs. (1) and (2) can be easily solved 
numerically, in particular in the co-propagating configuration and undepleted pump regime, which 
is of interest to us here. In this case the z-dependent elements on the diagonals of the SPM, XPM 
and Raman matrices, Jg, Jx, and J r, are obtained as previously discussed. 

When doing this, we found that both SPM and XPM effects have virtually no impact on the 
performance of Raman polarizers operating in the undepleted pump regime. In contrast, the form 
of the Raman matrix is of paramount importance. The larger the coefficients on the diagonal, the 
stronger the PDG. For moderate values of the polarization mode dispersion (PMD) coefficient, 
Raman diagonal terms only take appreciable values near the fiber input, as illustrated in Fig. 1. 
Therefore, the power of the pump beam is to be high, in order to provide significant amplification 
over the first few hundreds meters of the fiber. 

For analyzing the performance of Raman polarizers we identify three characteristic quantities: 
the degree of polarization (DOP) of the outcoming signal beam, its SOP, and the overall signal 
gain. The DOP and SOP characteristics are illustrated in Fig. 2. Since the signal SOP depends on 
the pump SOP, it is reasonable to define a quantity that measures the relative difference between 
these two SOPs. As usual, such quantity is the alignment parameter 

A K = ^ra ' (3) 

which is the cosine of the angle between the pump and the signal Stokes vectors, averaged over the 
ensemble of beams with random SOPs which models the unpolarized signal beam. The hypothesis 
that the signal SOP is attracted to the pump SOP is rooted in the model of isotropic fibers, in 
which Jjii = Jr2 = Jr2 = 1- In randomly birefringent fibers, the equality and even positivity of 
the three elements is not always the case, as exemplified in the plot of Fig. 1. In these cases, it is 
remarkable that the signal SOP is attracted to an SOP which is different from that of the pump. In 
spite of this observation, we found that for ideal Raman polarizers (those with DOP> 0.9), and in 
the range of lengths 0.001 < Lb < 0.05 and 0.0001 < L c < 0.05, given here in km, the signal SOP 
on average is attracted to the pump SOP, see Fig. 2. This is not the case in the counter-propagating 
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configuration, for which the appropriate alignment parameter Ah is different from that given in 
Eq. (3), see [5]. Moreover, the performance of Raman polarizers (namely, DOP) sensitively depends 
on the pump SOP, as demonstrated in Fig. 3. 




Fig. 1. Elements of the Raman matrix (Jri - black solid, Jfq - red dashed, and 
Jr3 ~ green dotted) as function of distance in the fiber for Lb(ui p ) = 0.016 km 
and L c = 0.05 km. (note that the black solid and red dashed curves coincide, i.e. 
Jri = Jr2-) 



Another important practical issue is the selection of fibers for Raman polarizers. The main 
parameter in this selection is the value of the PMD coefficient. In this respect we found that 
for obtaining a signal DOP close to unity (i.e., > 0.99) the PMD coefficient should be less than 
0.0145 ps/Vkm for, say, 8 W of pump power (as in Ref. [1]). Nevertheless we found that the PMD 
coefficient does not always provide full information about the fiber. For example in Fig. 2(d) we 
can see that two fibers with equal PMD coefficients exhibit a different performance as Raman 
polarizers. In one case, the DOP is 0.25, in the other - 0.45. For this reason, it is preferable to 
consider the beat and correlation lengths separately, rather than combining them into the single 
PMD coefficient, which for our model is expressed as [4]: D p = 2^/2n^fl7 c / \Lboj s ). 

The third characteristic of Raman polarizers is Raman gain, see Fig. 4. Even for a 1.5 km long 
fiber with 8 W of pump power we may have an enormous 55 dB gain that is almost twice the gain of 
the same Raman amplifier, but with a high value of the PMD coefficient. This means that Raman 
polarizers are simultaneously very efficient Raman amplifiers. Such values of gain are obtained in 
the undepleted regime, i.e. for input signal powers in the /iW range. For the mW range which is 
typical of telecom applications the analysis necessarily enters the depleted pump regime, to which 
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Fig. 2. DOP of the signal beam (black, solid) and alignment parameter Ay^ (red, 
dashed) as function of correlation length L c for the four SOPs of the pump beam: a) 
(l/\/3)(l, 1, 1); b) (1, 0, 0); c) (0, 1, 0); d) (0, 0, 1). Here and in Figs. 3 and 4, the 
value of the beat length Lb(oj p ) is indicated on the plots in km. The two ellipses on 
plot d) indicate one (of infinitely many) pair of points with equal PMD coefficients. 
Other parameters are (also used in Figs. 3, 4): input signal power 1 /mi; input pump 
power 8 W; g = 0.6 (W-krn)" 1 ; 7 = 1 (W-krn)" 1 ; a = 0.2 dB/km; L = 1.5 km. 
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Fig. 3. DOP of the signal beam for two SOPs of the pump beam which either 
maximize (black, solid) or minimize (red, dashed) the signal DOP. For each value 
of L c we perform a separate search for these two SOPs. 
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Fig. 4. Average Raman polarizer gain as a function of the correlation length. The 
pump SOP is (1, 0, 0) and the signal beam is initially unpolarized. 
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our theory can also be readily applied. 

In conclusion, we presented a theory for describing the interaction of two optical beams in ran- 
domly birefringent fibers via Kerr and Raman effects, and applied it to the quantification of the 
performance of Raman polarizers. 

We thank L. Palmieri for valuable comments. This work was carried out in the framework of the 
"Scientific Research Project of Relevant National Interest" (PRIN 2008) entitled "Nonlinear cross- 
polarization interactions in photonic devices and systems" (POLARIZON), and in the framework 
of the 2009 Italy-Spain integrated action "Nonlinear Optical Systems and Devices" (HI2008-0075). 
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